Can GJ 876 host four planets in resonance?

Prior to the detection of its outermost Uranus-mass object, it had been suggested that GJ 876 could host an Earth-sized planet in a 15-day orbit. Observation, however, did not support this idea, but instead revealed evidence for the existence of a larger body in a ~125-day orbit, near a three-body resonance with the two giant planets of this system. In this paper, we present a detailed analysis of the dynamics of the four-planet system of GJ 876, and examine the possibility of the existence of other planetary objects interior to its outermost body. We have developed a numerical scheme that enables us to search the orbital parameter-space very effectively and, in a short time, identify regions where an object may be stable. We present details of this integration method and discuss its application to the GJ 876 four-planet system. The results of our initial analysis suggested possible stable orbits at regions exterior to the orbit of the outermost planet and also indicated that an island of stability may exist in and around the 15-day orbit. However, examining the long-term stability of an object in that region by direct integration revealed that the 15-day orbit becomes unstable and that the system of GJ 876 is most likely dynamically full. We present the results of our study and discuss their implications for the formation and final orbital architecture of this system.     

A Minimal Spanning Tree algorithm for source detection in γ-ray images

We developed a source detection algorithm based on the Minimal Spanning Tree (MST), that is a graph-theoretical method useful for finding clusters in a given set of points. This algorithm is applied to γ-ray bi-dimensional images where the points correspond to the arrival direction of photons, and the possible sources are associated with the regions where they clusterize. Some filters to select these clusters and to reduce the spurious detections are introduced. An empirical study of the statistical properties of MST on random fields is carried out in order to derive some criteria to estimate the best filter values. We also introduce two parameters useful to verify the goodness of candidate sources. To show how the MST algorithm works in practice, we present an application to an EGRET observation of the Virgo field, at high Galactic latitude and with a low and rather uniform background, in which several sources are detected.     

SESTET: A spatially explicit stream temperature model based on equilibrium temperature

Stream-water temperature is a key variable controlling chemical, biological, and ecological processes in freshwater environments. Most models focus on a single river cross-section; however, temperature gradients along stretches and tributaries of a river network are crucial to assess ecohydrological features such as aquatic species suitability, growth and feeding rates, or disease transmission. We propose SESTET, a deterministic, spatially explicit stream temperature model for a whole river network, based on water and energy budgets at a reach scale and requiring only commonly available spatially distributed datasets, such as morphology and air temperature, as input. Heat exchange processes at the air–water interface are modelled via the widely used equilibrium temperature concept, whereas the effects of network structure are accounted for through advective heat fluxes. A case study was conducted on the prealpine Wigger river (Switzerland), where water temperatures have been measured in the period 2014–2018 at 11 spatially distributed locations. The results show the advantages of accounting for water and energy budgets at the reach scale for the entire river network, compared with simpler, lumped formulations. Because our approach fundamentally relies on spatially distributed air temperature fields, adequate spatial interpolation techniques that account for the effects of both elevation and thermal inversion in air temperature are key to a successful application of the model. SESTET allows the assessment of the magnitude of the various components of the heat budget at the reach scale and the derivation of reliable estimates of spatial gradients of mean daily stream temperatures for the whole catchment based on a limited number of conveniently located (viz., spanning the largest possible elevation range) measuring stations. Moreover, accounting for mixing processes and advective fluxes through the river network allows one to trust regionalized values of the parameters controlling the relationship between equilibrium and air temperature, a key feature to generalize the model to data-scarce catchments.     

Geochemical evidence of aeolian deposits in European soils

Hafnium (Hf) and zirconium (Zr) concentrations measured in over 4100 agricultural soil samples from Europe were assessed with the focus on their relationship to the distribution of aeolian deposits, such as loess and coversands. Comparison of extractable (aqua regia; ICP‐MS) and total (XRFS) concentrations shows that only 1.0 to 1.7% of the total Hf and Zr is chemically extractable because of the resistant nature of their host minerals. Resistate minerals, such as zircon, are commonly found in the predominantly silty fraction of loess deposits. In this study a statistical analysis of total Hf and Zr soil data from areas with and without loess was carried out to derive threshold values of 10 mg kg^?1 (Hf) and 318 mg kg^?1 (Zr). These values were subsequently applied across the project area in an attempt to indicate the presence of aeolian deposits. The spatial distribution of above‐threshold concentrations suggests a more extensive and coherent loess belt across central and eastern Europe, providing additional evidence of loess across Brittany, Aquitaine and near the Vosges mountains in France as well as in the basins of central and northern Spain. Above‐threshold concentrations were also detected in loess regardless of its thickness, emphasising the importance of the abundance of zircon in the upper part of the soil profile rather than the actual thickness of the deposit. Soil data however, failed to indicate various loess facies within most of the Pannonian Basin, suggesting that this approach only works where deposits contain a sufficient amount of zircon. This may also explain why the extensive coversands across northern Germany and Poland were largely undetected. This study demonstrates that continental‐scale soil geochemical data can help identify and map the distribution of zircon‐rich loess and coversand, and subsequently enhance and improve current knowledge of the extent of these deposits.     

Nonstationary stochastic response of structural systems equipped with nonlinear viscous dampers under seismic excitation

Nonlinear viscous dampers are supplemental devices widely used for enhancing the performance of structural systems exposed to seismic hazard. A rigorous evaluation of the effect of these damping devices on the seismic performance of a structural system should be based on a probabilistic approach and take into account the evolutionary characteristics of the earthquake input and of the corresponding system response. In this paper, an approximate analytical technique is proposed for studying the nonstationary stochastic response characteristics of hysteretic single degree of freedom systems equipped with viscous dampers subjected to a fully nonstationary random process representing the seismic input. In this regard, a stochastic averaging/linearization technique is utilized to cast the original nonlinear stochastic differential equation of motion into a simple first‐order nonlinear ordinary differential equation for the nonstationary system response variance. In comparison with standard linearization schemes, the herein proposed technique has the significant advantage that it allows to handle realistic seismic excitations with time‐varying frequency content. Further, it allows deriving a formula for determining the nonlinear system response evolutionary power spectrum. By this way, ‘moving resonance’ effects, related to both the evolutionary seismic excitation and the nonlinear system behavior, can be observed and quantified. Several applications involving various system and input properties are included. Furthermore, various response parameters of interest for the seismic performance assessment are considered as well. Comparisons with pertinent Monte Carlo simulations demonstrate the reliability of the proposed technique. Copyright © 2014 John Wiley & Sons, Ltd.     

Response reduction factor of irregular RC buildings in Kathmandu valley

Most current seismic design includes the nonlinear response of a structure through a response reduction factor(R). This allows the designer to use a linear elastic force-based approach while accounting for nonlinear behavior and deformation limits. In fact, the response reduction factor is used in modern seismic codes to scale down the elastic response of a structure. This study focuses on estimating the actual ‘R' value for engineered design/construction of reinforced concrete(RC) buildings in Kathmandu valley. The ductility and overstrength of representative RC buildings in Kathmandu are investigated. Nonlinear pushover analysis was performed on structural models in order to evaluate the seismic performance of buildings. Twelve representative engineered irregular buildings with a variety of characteristics located in the Kathmandu valley were selected and studied. Furthermore, the effects of overstrength on the ductility factor, beam column capacity ratio on the building ductility, and load path on the response reduction factor, are examined. Finally, the results are further analyzed and compared with different structural parameters of the buildings.     

Travertine deposition and faulting: the fault-related travertine fissure-ridge at Terme S. Giovanni,Rapolano Terme (Italy)

In this paper we describe an example of travertine fissure-ridge development along the trace of a normal fault with metre displacement, located in the eastern margin of the Neogene–Quaternary Siena Basin, in the Terme S. Giovanni area (Rapolano Terme, Italy). This morphotectonic feature, 250 m long, 30 m wide and 10 m high, formed from supersaturated hot waters (39.9°C) flowing from thermal springs aligned along the trace of the normal fault dissecting travertines not older than Late Pleistocene (24 ± 3 ka). A straight, continuous fissure with a maximum width of 20 cm occurs at the top of the ridge, along its crest. Hot fluids flow from cones mainly located at the extremities of the ridge, where travertine is depositing. The travertine fissure-ridge shows an asymmetrical profile since it buries the fault scarp. The difference in height of slopes corresponds to the vertical displacement of the normal fault. Fissuring of the recent travertine deposits along the strike of the crestal fissure, as well as recent hydrothermal circulation, lead us to believe that the Terme S. Giovanni normal fault may be currently active. On the whole, the Terme S. Giovanni fissure ridge can be defined as a travertine fault trace fissure-ridge, adding a helpful example for studying the relationship between faulting and travertine deposition.     

Seismic driven probabilistic classification of reservoir facies for static reservoir modelling: a case history in the Barents Sea

In this paper we present a case history of seismic reservoir characterization where we estimate the probability of facies from seismic data and simulate a set of reservoir models honouring seismically‐derived probabilistic information. In appraisal and development phases, seismic data have a key role in reservoir characterization and static reservoir modelling, as in most of the cases seismic data are the only information available far away from the wells. However seismic data do not provide any direct measurements of reservoir properties, which have then to be estimated as a solution of a joint inverse problem. For this reason, we show the application of a complete workflow for static reservoir modelling where seismic data are integrated to derive probability volumes of facies and reservoir properties to condition reservoir geostatistical simulations. The studied case is a clastic reservoir in the Barents Sea, where a complete data set of well logs from five wells and a set of partial‐stacked seismic data are available. The multi‐property workflow is based on seismic inversion, petrophysics and rock physics modelling. In particular, log‐facies are defined on the basis of sedimentological information, petrophysical properties and also their elastic response. The link between petrophysical and elastic attributes is preserved by introducing a rock‐physics model in the inversion methodology. Finally, the uncertainty in the reservoir model is represented by multiple geostatistical realizations. The main result of this workflow is a set of facies realizations and associated rock properties that honour, within a fixed tolerance, seismic and well log data and assess the uncertainty associated with reservoir modelling.